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Abstract. We construct discrete analogues of Dirac structures by considering the geometry of symplectic 
maps and their associated generating functions, in a manner analogous to the construction of continuous 
Dirac structures in terms of the geometry of symplectic vector fields and their associated Hamiltonians. We 
demonstrate that this framework provides a means of deriving implicit discrete Lagrangian and Hamilton- 
ian systems, while incorporating discrete Dirac constraints. In particular, this yields implicit nonholonomic 
Lagrangian and Hamiltonian integrators. We also introduce a discrete Hamilton-Pontryagin variational 
principle on the discrete Pontryagin bundle, which provides an alternative derivation of the same set of in- 
tegration algorithms. In so doing, we explicitly characterize the discrete Dirac structures that are preserved 
by Hamilton-Pontryagin integrators. In addition to providing a unified treatment of discrete Lagrangian 
and Hamiltonian mechanics in the more general setting of Dirac mechanics, it provides a generalization of 
symplectic and Poisson integrators to the broader category of Dirac integrators. Since discrete Lagrangians 
and discrete Hamiltonians are essentially generating functions of different types, the theoretical framework 
described in this paper is sufficiently general to encompass all possible Dirac integrators through an appro- 
priate choice of generating functions. 
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1. Introduction 

Dirac structures, which can be viewed as simultaneous generalizations of symplectic and Poisson struc- 
tures, were introduced in Courant In the context of geometric mechanics [H; [H; EH > Dirac structures 
are of interest as they can directly incorporate Dirac constraints that arise in degenerate Lagrangian sys- 
tems [lij , interconnected systems [8[ , and nonholonomic systems Q , and thereby provide a unified geometric 
framework for studying such problems. 

From the Hamiltonian perspective, these systems arc described by implicit Hamiltonian systems, and a 



comprehensive review of Dirac structures in this setting can be found in Dalsmo and van der Schaft 11|. 
This approach is motivated by earlier work on almost-Poisson structures that describe nonholonomic systems 
using brackets that fail to satisfy the Jacobi identity j(| [27| ■ In the context of systems with symmetry, Dirac 
analogues of symplectic Q and Poisson [13] reduction have been developed. 



l 



On the Lagrangian side, degenerate, interconnected, and nonholonomic systems can be described by 
implicit Lagrangian systems, which were introduced in the context of Dirac structures in Yoshimura and 
Marsden (2q |. The corresponding variational description of implicit Lagrangian systems was developed in 
[2^], with the introduction of the Hamilton-Pontryagin principle on the Pontryagin bundle TQ®T*Q, which 
yields the Legcndrc transformation, as well as Hamilton's principle for Lagrangian systems and Hamilton's 
phase space principle for Hamiltonian systems. Implicit Lagrangian systems with constraints and external 
forces are described by the Lagrange-d'Alcmbcrt-Pontryagin principle, and in the presence of symmetries, 
the corresponding Lagrangian reduction theory was developed in [30(, and the cotangent bundle reduction 
theory in [3l| . Generalizations of Dirac manifolds to Dirac anchored vector bundles provide a category that 
is closed under Dirac reduction, which provides the necessary setting for Dirac reduction by stages 



In the context of geometric numerical integration |14t |2l| , which is concerned with the development of 
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numerical methods that preserve geometric properties of the corresponding continuous flow, variational in- 
tegrators that preserve the symplectic structure can be systematically derived from a discrete Hamilton's 



principle [251 ]. and can be extended to asynchronous variational integrators 23j that preserve the multi- 
symplectic structure of Hamiltonian partial differential equations. The discrete variational formulation of 
Hamiltonian mechanics was developed in [2(| as the dual, in the sense of optimization, to discrete Lagrangian 
mechanics. Discrete analogues of the Hamilton-Pontryagin principle were introduced in [B|; [l8[ for particular 
choices of discrete Lagrangians. 

Contributions of this paper. In this paper, we introduce discrete analogues of Dirac structures, and 
show how they describe discrete implicit Lagrangian and Hamiltonian systems. The construction relies on 
the observation that the continuous Dirac structures used to define implicit Lagrangian systems arise from 
geometric properties of infmitesimally symplectic vector fields. By analogy, we construct discrete Dirac struc- 
tures that are derived from properties of symplectic maps, and demonstrate that they yield implicit discrete 
Lagrangian and Hamiltonian systems, and recover nonholonomic integrators that are typically derived from 
a discrete Lagrange-d'Alcmbcrt principle. 

We also introduce a discrete Hamilton-Pontryagin principle on the discrete Pontryagin bundle (Q x Q) © 
T*Q, that provides a variational characterization of implicit discrete Lagrangian and Hamiltonian systems 
that we previously described in terms of discrete Dirac structures, and which reduce to the standard varia- 
tional Lagrangian integrators [25| and variational Hamiltonian integrators [2(|. Furthermore, we introduce 
a discrete Lagrange-d'Alcmbcrt-Pontryagin principle to allow the incorporation of Dirac constraints, and 
which recover nonholonomic integrators. We also describe a discrete Hamilton's phase space principle, which 
provides a variational formulation of discrete Hamiltonian mechanics, and a discrete Hamilton-d'Alcmbcrt 
principle in phase space, which addresses the issue of constraints. 

In addition to providing a characterization of the discrete geometric structure that is preserved by 
Hamilton-Pontryagin integrators, we also characterize the corresponding discrete variational principles in 
an intrinsic manner that provides a theory of discrete Dirac mechanics that is valid semi-globally on the 
discrete Pontryagin bundle, i.e., on the preimage of a neighborhood of the diagonal of Q x Q. We also 
provide a correspondence between discrete Lagrangian and Hamiltonian mechanics by introducing a discrete 
generalized Legcndrc transformation that is valid even without the assumption of hypcrrcgularity. 

Outline of this paper. The paper is organized as follows. The first part of the paper is concerned with 
the geometry of discrete Dirac mechanics. In Section 2, we review the theory of generating functions of 
symplectic maps, and derive discrete maps that we will use in the construction of discrete Dirac structures. 
In Section 3, we describe the corresponding theory of infmitesimally symplectic vector fields, and derive 
continuous maps that describe continuous Dirac structures. In Section 4, we review the continuous theory of 
Dirac structures and implicit Lagrangian and Hamiltonian systems, and construct a corresponding discrete 
theory. 

The second part of the paper addresses the discrete variational structure of discrete Dirac mechanics. 
In Section 5, we review the Hamilton-Pontryagin principle and implicit Lagrangian systems. In Section 6, 
we introduce the discrete generalized Legcndrc transformation. In Section 7, we introduce both the local 
and intrinsic discrete Hamilton-Pontryagin principle, and show that they yield implicit discrete Lagrangian 
and Hamiltonian systems. In Section 8, we introduce the local and intrinsic discrete Lagrange-d'Alembert- 
Pontryagin principle, which incorporates discrete constraints that can model both interconnections and 
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nonholonomic constraints. In Sections 9 and 10, we consider the corresponding local and intrinsic discrete 
variational principles on the Hamiltonian side, which are respectively the discrete Hamilton's principle on 
phase space, and the discrete Hamilton-d'Alembcrt principle on phase space. In Section 11, we provide some 
concluding remarks and future directions. 

2. The Geometry of Generating Functions 

2.1. Generating Functions. Let us first review the theory of generating functions following [l[ and 24|. 
The key fact is the following (Proposition 5.2.1 of [l[): 

Proposition 2.1. Let (Pq, fig) and (Pi, fii) be symplectic manifolds, and ni : Pq x Pi — > Pj be the projections 
onto Pi for i = 0,1, and let us define f2p oX Pi £ A (-^0 x Pi) as follows: 

^PoxPi := 7r*fii — 7rQf2o. (2.1) 

Then 

(i) f2p x p x is a symplectic form on Pq x Pi. 

(ii) A map F : Pq — > Pi is symplectic if and only if ipQp a *p 1 = 0, where ip : Tp — > Pq x Pi is an 
inclusion and Tp is the graph of F. 

In particular we are interested in the case (Pi,Oj) = (T*Qi,VLi) with the canonical symplectic form 
n t = -d@i on T*Qi for i = 0, 1. Then 

^•Qoxt-q, := tt^! - n*n = <(-d9i) - ^(-d8o) = -d^Gi - ^6 ) = -de^oxr^. (2.2) 

where we defined Q T 'Q xt*q 1 £ A 1 ( T *Qo X ^*Qi) by 

©T*OoxT«Oi : = 7r*0i - 7I"o©0- (2.3) 

So 

ipOT»Q oX T*Qi = *p(-rf@T*QoxT*Q!) = -^(«F@T«Q xT'Qi )• (2.4) 

Therefore Part Jip of the above proposition now reads 

F : T*Q Q -> T*Q 1 is symplectic ^pe^QoxT-Qj - 0. (2.5) 

However, by the Poincare lemma, locally d(«p©T*Q xT*Qi) = if and only if ipQr^QayT'Qx = dS for some 
function S : Tp — > R. Such a function S* is called a generating function. To summarize, 

F : T*Q -> T*Qj is symplectic '^. y ip6 T *Q xt-q 1 = dS for some function S : T F -> R. (2.6) 

2.2. Parametrization of Tp. Suppose there exists a manifold M and a map ip^ : M — > Tp that gives a 
(local) parametrization of the graph Tp. Now define S : M — ► R to be S := <p* M S = S o <p^/. Then we can 
restate Eq. (|2.6[) in the following way: 

F : T*Q -> T*Qi is symplectic (i^ )*6 t .q oXT .q 1 = for some function S : M -> R, (2.7) 

where i£f : A/ — > T*<2o x T*Qi is defined by := i F o 

For the remainder of this section, we consider the special case with Qq = Qi = Q. This is a natural 
setting for doing mechanics since in this case F : T*Q — > T*Q describes a flow on the cotangent bundle T*Q. 
We choose three different parametrization based on the classification given by Goldstein et al. 13 1. 

2.3. Generating Function of Type 1 and the Map 7^ : T*Q x T*Q -^T*(Q x Q). Let M be Q x Q. 

Then the flow F on PQ is symplectic if and only if there exists Si : Q x Q — > R such that 

^ XQ )*e T .QxT.Q=dSi. (2.8) 

Suppose F : (<Zo,J>o) l— ^ (qi , _Pi ) , or cquivalcntly, 

^ xQ :Qxg^rQxrg ; ( qo , qi ) ^ (qupi)), (2-9) 

where pq and pi are considered to be functions of go and qy. Then we can write Eq. (|2.8[) as follows: 

Pidqi -p dq = DySydqo + D 2 Sid qi , (2-10) 

which gives 

po = -DiSi, pi = D 2 Si. (2.11) 
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This gives rise to a map Kq : T*Q x T*Q — ► T*(Q x Q) so that the diagram 

T*Q x T*Q ^ T*(Q x Q) 



,QxQ 



Q x Q 



commutes. To be more specific, one has 

((go,Po),(gi,Pi)) 



(qo,qi,DiSi,D 2 Si) 



(2.12a) 



(<?o,<7i) 



In view of Eq. (|2.1ip . we obtain 



K o : ((go,po), (gi, pi)) !->• (go, gi, -po,pi)- 



(2.12b) 



(2.13) 



2.4. Generating Function of Type 2 and the Map il b d+ : T*Q x T*Q -> Let M be H+, whose 

local coordinates are (go,Pij]l Then the flow F on T*Q is symplectic if and only if there exists 52 : TC+ — > K 
such that 



(v + )*e^ 



T'QxT* 



where 



e (2) 

T*QxT*Q 



= d(qipi) - Ot*QxT*q = Podqa + q\dp x . 



(2.14) 
(2.15) 



( 2) 

Note that using Qt,q x t*q m place of Ot*qxt*q does not affect the argument outlined in Section |2~21 since 

d®T*QxT*Q = -dO T *QxT*Q = &<T*QxT*Q- 

Suppose F : (go,Po) l— * (gi,Pi), or equivalcntly, 



V + : -> T*g x T*Q; (g ,Pi) ^ ((go,Po), (gi,Pi)), 



where po an d 9i are considered to be functions of q$ and p\. Then we can write Eq. (|2. 14)1 as follows: 

p dq + qidpi = D 1 S 2 dq + D 2 S 2 dp 1 , 

which gives 

Po = D 1 S 2 , qx = D 2 S 2 . 
This gives rise to a map Q b d+ : T*Q x T*Q — > T*H+ so that the diagram 



(2.16) 

(2.17) 
(2.18) 



T*Q x T*Q 



■ T*n+ 



dS 2 



(2.19a) 



commutes. To be more specific, one has 

((go.Po), (gi,Pi)) 



(go, Pi, DiS 2l D 2 S 2 ) 



In view of Eq. (|2.18p . we obtain 



(go, Pi) 



: ((go,po), (gi,Pi)) >->• (go,pi,po,gi) 



(2.19b) 



(2.20) 



■*"We can think of 7i+ as a submanifold of Q X T*Q with local coordinates (qo, (gi,pi)) where 51 is dependent on qo and pi. 
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2.5. Generating Function of Type 3 and the Map fl h d _ : T*Q x T*Q -> T*H-. Let M be W_, whose 
local coordinates are {pi, goU- Then the flow F on T*Q is symplectic if and only if there exists S3 : — > R 
such that 

(%-ye% xT . Q = ds 3 , (2.21) 

where 

®t*qxt*q '■= - d (qoPo) - Qt'Qxt*q = -qodpo - Pidqi. (2.22) 

Again, we can use &t* QxT ,q in place of &t'Qxt*q since cIQt} QxT , q = -d<d T * QxT 'Q = ^t*qxt*q- 
Suppose F : (qo,Po) (<7i, Pi), or equivalently, 

: H_ -> T*Q x T*Q; (po,?i) 1- ((fflj.Do), (?i,Pi)), (2.23) 

where go an d Pi are considered to be functions of po an d q\. Then wc can write Eq. (|2.21[) as follows: 

- q dp - pidqi = DiSzdp + D 2 S z dqi 1 (2.24) 

which gives 

q = Pl = -D 2 S 3 . (2.25) 
This gives rise to a map f2^_ : T*Q x T*Q — > T*H- so that the diagram 



T*Q x T*Q 



■T*U- 



dS 3 



H. 



commutes. To be more specific, one has 

{(Qo,Po), (<7i,Pi)) 



(po,qi,D 1 S 3l D 2 S 3 ) 



(2.26a) 



(2.26b) 



In view of Eq. (|2.25[) . wc obtain 



(Po,3i) 



((9o,po),(9i,pi)) !-»• (po, qi, -qo,-px)- 



(2.27) 



3. Symplectic Flows as the Infinitesimal Limit of Symplectic Maps 

We can "infinitesimalize" the above discussions to recover the familiar notions of symplectic flows, Hamil- 
tonian and Lagrangian systems, and also the maps kq : TT*Q — » T*TQ and tt b : TT*Q — » T*T*Q that are 
found in Yoshimura and Marsden [2^ . 

3.1. Hamiltonian Flows and the Map fi b : TT*Q -> T*T*Q. The key idea is to regard the flow 
Fx : T*Q — > T*Q of a vector field X e X(T*Q) as the infinitesimal limit of the above discussions of 
symplectic maps. In fact the following definition of symplecticity of the flow Fx is analogous to Eq. (|2.5|l : 



Definition 3.1. The flow F x : T*Q -> T*Q of a vector field X G X{T*Q) is called symplectic if = 0. 

By Cartan's magic formula, and since f2 is a closed two-form, i.e., dfl = 0, a symplectic vector field X 
can be equivalently characterized by the property, 

= £ x tt = ix(dCt) + d(i x Q) = d(i x ty, (3.1) 

which is to say that a vector field X S X(T*Q) is symplectic if i x Q, is closed, i.e., d(i x Q) = 0. 
Now again by the Poincare lemma, locally we can restate this as follows: 

F x : T*Q -> T*Q is symplectic '^- y i x fl = dH for some function H : T*Q -> R, (3.2) 



We can think of H— as a submanifold of T*Q X Q with local coordinates ((<JOiPo)i 91 ) where go is dependent on po and q\. 
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which is again analogous to Eq. (|2.6[) . So the Hamiltonian H : T*Q — ► R is an infinitesimal analogue of 
generating functions. Furthermore, the above local statement leads to the well-known global definition of the 
Hamiltonian flows: 

Definition 3.2. The flow Fx : T*Q -> T*Q of a vector field X e X(T*Q) is called Hamiltonian if is 
exact, i.e., ix H = d# for some function H : T* Q ^ M.. 

Now analogously to Section [21 we can define $l b so that 

52° 



TT*Q 



.T*T*( 



X 



dH 



(3.3a) 



commutes. To be more specific, one has 



(q,p,dH/dq,dH/dp) 



(3.3b) 



{q,p) 

Note that, in coordinates, we can write ix^l = dH as follows: 

dH . dH 
dp 1 P dq ' 
So in view of this set of equations, we obtain 

ft b : (q,p, Sq, 5p) h-> (q,p, ~8p, Sq). 



q 



(3.4) 



(3.5) 



3.2. Lagrangian Flows and the Map kq : TT*Q — > T*TQ. We consider the Lagrangian analogue of the 
above construction. Given the Legendre transformation FX : TQ — > T*Q, we construct Ql = FX*f2, and 
consider a second-order vector field Xl S X(TQ) that preserves the Lagrangian symplectic form, i.e., 

£ XL toL = 0. (3.6) 

We introduce the Lagrange one-form 9l = FX*6 = Q^jdq. Then, Ql = ~ rf@L, and 

= £ Xl Vl = £x L (-dO L ) = -d£x L Q L , (3.7) 

as the Lie derivative commutes with the exterior derivative. Since £x L is closed, by the Poincare lemma, 
there exists a local function L : TQ — > M, such that, 

£x L ®L=dL. (3.8) 

This is the intrinsic Euler-Lagrange equation expressed in terms of the Lagrangian (Section 3.4.2 of [3]), 
and is equivalent to the intrinsic Euler-Lagrange equation written in terms of the energy function (Equation 
(7.3.5) of 0), 

i XL n L =dE 1 (3.9) 
as the following discussion demonstrates. By applying Cartan's magic formula, we obtain 

dL = £ Xl 0l =ix L {dQ L ) + d{ix L @L) = -%x^l + d {ix L (FX*6)) , (3.10) 

which implies 

ix L Vi L = d{ix L {n*@)-L). (3.11) 

The expression in the parentheses on the right hand side is precisely the intrinsic expression for the energy 
function E : TQ — > R, as the following coordinate computation shows: 

dL I dh \ 

E(q, v) = —v- L(q, v) = ( -^dq, (q, v) ) - L(q, v) = (i Xl (¥L*Q) - L) (q, v), (3.12) 
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where we used the fact that Xl is second-order, i.e., q = v. As such, (|3.8[) is equivalent to (|3.9[) . In 
coordinates, we can write £x L ®L — dL, expressed in terms of the FL-related vector field X G X(T*Q), i.e., 
Xo¥L = TFL oX L ,as follows: 



Now we can define kq so that 



P 



dL 

dv ' 



v, 



dL 



(3.13) 




T*Q 

commutes. To be more specific, one has 



FL 



TQ 



T*TQ 




(3.14a) 



In view of Eq. (|3.13[) . we obtain 



(q, v, q, v) (q,v,dL/dq,dL/dv) 



(3.14b) 



K Q ■ (q,P,Sq,Sp) i-> (q,6q,5p,p). 



(3.15) 



4. Discrete Dirac Structures and Implicit Discrete Lagrangian Systems 

The maps Kq and f2^ ± defined in Eqs. (|2 . 13[> . (|2.20p . and (|2.27[) provide a discrete counterpart of the 
framework for (continuous) Dirac mechanics developed by Yoshimura and Marsden 2$k 22| ■ Considering the 
fact that the discrete Lagrangian and Hamiltonians are generating functions of Types 1, 2, and 3 [20| . this 
is a natural setting for the implicit discrete Lagrangian and Hamiltonian systems, as we shall see. 



For the purpose of comparison, let us first briefly review continuous Dirac mechanics following [281 : |2 



4.1. Continuous Dirac Mechanics. 
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(q,p,-Sp,8q) 



{q,p,6q,5p) 



(q, Sq,6p,p) 




(4.1b) 



{q,Sq) © (q,p) 



{q,Sq) 



4.1.2. Symplectic Forms. The induced symplectic one-forms on TT*Q are 

x ■■= (n b ye T , T * Q = -s P dq + s q dp 7 x:=(K Q ye T -TQ = Spdq+pd{Sq). (4.2) 

and 

fl TT , Q := -d\ = d\ = dqA d(6p) + d(Sq) A dp. (4.3) 

4.1.3. Dirac Structure and Implicit Lagrangian System. A distribution Aq C TQ induces a Dirac structure 
on T*Q as follows: 

D Aq {z) := {(v z ,a z ) G T Z T*Q x T 2 *T*Q | v z G A T . Q (z), a*K) = fl(v z ,w z ) for u> 2 G A T , Q (z)} , (4.4) 
or equivalently, 

£>A Q (z) := {(^,«.) € T Z T*Q x T 2 *T*Q | i; 2 G A T . Q (z), a,-ft k («„) G A^, Q (z)} (4.5) 

where A T . Q := (Tttq^^Aq). 

Let 7 Q := fi b o (kq)" 1 : T*TQ -> T*T*Q. For a given Lagrangian L : TQ -> R, define DT := 7q o dT. 
Let X G X(T*Q), i.e., X : T*Q — > TT*Q. Then an implicit Lagrangian system (L,Aq,X) is defined as 
follows: 



(X,DL) G D Ac 



This gives the implicit Lagrangian system: 



* / \ dL dL . . , . 

g = veA (g), p=— , p -_eA£(?). 



(4.6) 
(4.7) 



Consider the special case Aq = TQ. Then 

Da q (z) := {(« z ,a 2 ) GT 2 T*Q xT 2 *T*Q | a z =Q b (v z )} 

Then the implicit Lagrangian system is written, in coordinates, as follows: 

dL dL 

q = v, p = -j-, p = 

aq ov 

4.1.4. Dirac Structure and Implicit Hamiltonian System. Given a Hamiltonian H : T*Q — > R, an implicit 
Hamiltonian system (H, Aq,X) is defined as follows: 

{X,dH)GD AQ , 

which gives the implicit Hamiltonian system: 



(4.8) 
(4.9) 



9H . , , . dH . „ , . 

9 = ^ eAQ(<z) ' p+ a7 eA « ((?) - 



(4.10) 
(4.11) 



Again, for the special case Aq = TQ, the implicit Hamiltonian system becomes 

dH . dH 
dq ' 

which is the standard Hamiltonian system. 



1 = "a - ' P 
op 



(4.12) 



4.2. (+)-Discrete Dirac Mechanics. 
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{qo,pi,po,qi) 



(qo,qi,-po,pi) 




(4.13b) 



4.2.2. Symplectic Forms. We can reinterpret the maps Kq and Q b d+ in connection to the discussions by 
Yoshimura and Marsden [28[ in the following way. The maps Kq and f2^ + induce two symplectic one-forms 
on T*Q x T*Q. One is 

Xd+ ■= (^ d+ )*6T*-H+ = Pa dq + qi dp 1 = Qt, QxT , q , (4.14) 

and the other is 

•W : = { k q)*@t*(QxQ) = -Po dq + pi dqi = Qt*qxt*q- (4.15) 
As shown above, they are related as follows: 

^T'QxT'Q = ~d\ d+ = dxd+ = dqi A dpi - dq A dp Q . (4.16) 

This nicely parallels with the corresponding discussions of the maps kq and fl b in [2^ and Section 14.1.21 

4.2.3. Discrete constraint distributions . We first introduce the discrete constraint distribution^ which we 
denote Ag C Q x Q, which is a submanifold of Q x Q with the property that it contains the diagonal of 
Q x Q, i.e., 

A = {(<?, g) | qeQ}c A d Q . (4.17) 
This discrete constraint distribution induces a continuous constraint distribution Aq C TQ as follows. 
Consider smooth curves on Q, endowed with the equivalence relation ip ~ ip <==^ y(0) = tp(0), Dip(0) = 
Dip(0). Given a curve ip £ C°°((— e, e), Q) with the properties <p(0) = q, Dip(0) = v, we identify the 
equivalence class of curves [ip] with the tangent vector v q € T q Q. Now, we introduce the class of curves on 
Q that are compatible with the discrete constraint distribution, 

C A * := {p e C°°((-e, c), Q) | Vr e (0, e), (ip(-r), <p(0)), (^(0), V (r)) G A Q } . (4.18) 

Then, the continuous constraint distribution Aq is defined by the property, 

<p G C A d =>> M G A Q . (4.19) 



We will adopt the notational convention that a continuous constraint distribution on M is denoted by Am C TM, and a 
discrete constraint distribution on M is denoted by A M C M X M. 
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The annihilator of Aq is denoted by Aq C T*Q, and is defined, for each q G Q, as 

Aq(<?) := K e T^Q I \/v q G A Q , <a 9 ,« g > = 0} . (4.20) 

Finally, these induce the discrete constraint distribution, A^.g := (ttq x 7Tq) _1 (Aq) C T*Q x T*<3, and the 
annihilator distribution on the discrete Pontryagin bundle, A^ ± := (f2^j.) (Aq x Aq) C T*H±, which are 
explicitly given by 

A^ Q = {((q ,p ),(q 1 ,p 1 )) GT*Q xT*Q \ (q , qi ) G A Q } , (4.21) 
A w+ = {(<7 ,Pi,Po,<7i) G T*H+ | p G A° Q (qo), Pl G Aq( 9i )} , (4.22) 
A w _ = {bo, 9i, -9o, -Pi) G | po E Aq( 9 o),Pi G Aq( 9i )} . (4.23) 

4.2.4. (+)-Discrete Dirac Structure and Implicit Discrete Lagrangian System. Let us define the (+)-discrete 
Dirac structure by 

D d A + Q (z):={((z,z 1 ),a z+ )€({z}xT*Q)xT; + H + \ (z, z 1 ) G A£. Q , a z+ - fi> + ((z, z 1 )) G A° H+ } , (4.24) 

where if z = (e?,p) and z 1 = (f? 1 ,?? 1 ) then z + := (g,^ 1 ). 

Let 7q + := o (kq)' 1 : T*(Q x Q) -> T*H + , and for a given discrete Lagrangian L d : Q x Q — > K, 
define := ^ odL. Now let us write = ((g°,p°), G x T*Q, and define an implicit 

discrete Lagrangian system (L d , Aq, X d ) as follows: 

n d+ (X%) e A° H+ . (4.25) 
DiL d ,ql), (4.26) 

(4.27) 

So we obtain the set of equations, 

p° k+1 = D 2 L d (ql,ql) G A° Q {q\), p° k + DM&ql) G A Q (q° k ), q l k = q° k+1 , (q° k ,q° k+1 ) G Aq, 

(4.28) 

which we shall call the implicit discrete Euler- Lagrange equations. 

Consider the special case Aq = Q x Q, which implies Aq = TQ. Then 

:= {((^W) e ({*> >< T *Q) x I = n d+ ((^^))} > (4-29) 

Then the implicit discrete Euler-Lagrangc equations have the form 

p k+1 =D 2 L d (q k ,ql), pl = -D 1 L d (ql,ql), q\ = q° k+v (4.30) 

4.2.5. (+)-Discrete Dirac Structure and Implicit (+)- Discrete Hamiltonian System. Let H d+ : TL + — > R be 
given. Then 

dH d+ {ql,pl) = D l H d+ {ql,p l k )dql + D 2 H d+ {ql,p 1 k )dp 1 k G T*H + (4.31) 
An implicit (+) -discrete Hamiltonian system (H d+ , Aq, X d ) is defined as 

(X*, dff^,^)) G D d A + (q° k ,q° k+1 ) G Aq and dH d+ - Vt d+ {X k d ) G A^ + . (4.32) 

Now 

dH d+ {ql,q l k ) = {ql,p 1 k ,D l H d+ ,D 2 H d+ ), (4.33) 

and 

n d+ (X k d ) = {ql,pl +1 ,pl,q° k+1 ). (4.34) 

So we obtain the set of equations 

p Q k - D.HMqlpl) G A° Q (q° k ), q° k+1 = D 2 H d+ {ql,p\), p\ - p° k+1 G A Q (ql), (q° k ,q° k+1 ) G A d Q , 

(4.35) 

which we shall call the implicit (+)-discrete Hamilton's equations. 



(XlV + L d (q° k ,ql)) G D d + Q ^ (q° k ,q° k+1 ) G A d Q and S+L, - 

Now 

S+L d (^, q\) = ^{qlqlDyLa, D 2 L d ) = {q a k ,D 2 L d , - 

and 

n d+ {X k d ) = {ql,pl +1 ,pl,ql +l ). 
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If Aq = Q x Q, the discrete Dirac structure is given by (|4.29[) . and the implicit (+)-discrete Hamilton's 
equations have the form 

pi = D 1 H d+ (q k ,p 1 k ), ti+i = D 2 H d+ {qlp\), pl=p° k+1 . 
Note that this is essentially the (+)-discrete Hamilton's equations in Lall and West [2fj| . 
4.3. (— )-Discrete Dirac Mechanics. 



(4.36) 



4.3.1. The Big Diagram. The maps Kq and fl d _ defined in Eqs. (|2.13[) and (|2 . 27[) give rise to the following 
diagram. 



T*H- 




QxQ 



T*(Q x Q) 



(4.37a) 



(po,9i, ~9o, -Pi) 



((9o,Po),(gi,pi)) 



(po, qi) •* (go, 9i ) © (po, qi) 



(90, gi) 



(9o, 9i, ~Po,Pi) 




(4.37b) 



Xd- ■= W d -)*®T*T>Q = -Pi dqi - 90 dp = ©r.Q xT .Q, 
A d - := (kq)*Q T '(QxQ) = -Podq + Pidq x = Qt'Qxt-q, 



4.3.2. Symplectic Forms. As in the (+)-discrete case, we can reinterpret the maps Kq and as follows. 
We have 

(4.38) 
(4.39) 
(4.40) 



and 

and then 



n T , QxT *Q := -d\d- = dxd- = dqi A dp\ - dq A dp . 

4.3.3. (-)-Discrete Dirac Structure and Implicit Discrete Lagrangian System. Let us define the (-)-discrete 
Dirac structure by 

Dt Q (*) : = { ((*. A ) e ({4 x T*Q) x r 2 *_ H_ I (z, z 1 ) G A*, Q , a z _- Q d _ ((z, z 1 )) eA° | (4.41) 

where if z = (q,p) and z 1 = (g 1 ,?? 1 ) then z- := (g 1 

Let 7g~ := o (kq)" 1 : T*(Q xQ)^ T*H-, and for a given discrete Lagrangian L d : Q x Q — > M, 
define Q~L d := 7q~ o cLL. Now define an implicit discrete Lagrangian system (L d , Aq, X d ) as follows: 

(A^STL^,^,)) G £>1- (gg,gg +1 ) G A$ and - fi^_(*$) G A° n _. (4.42) 

Now 

JT^fag,^) =J Q ~(qlq 1 k ,D 1 L d ,D 2 L d ) = (-D 2 L d , q° k+1 , -q° k , -D x L d ), (4.43) 

and 

n b d _(X k d ) = (Plq k+1 ,- q ° k ,-p k+1 ). (4.44) 
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So we recover the implicit discrete Euler-Lagrange equations 

P° k+1 = D 2 L d (q° k ,ql) e A° Q (ql), p° + DM&Qk) e A^(gg), - q° k+1 , (q k ,q° k+1 ) S Aq, 

(4.45) 

that we previously obtained in (|4.28|) . 
If A Q = Q x Q, then 

Dt Q (z) := {((*, G ({4 x T*Q) x T^Ti- \ a z _ = fi>_ ((z.z 1 ))} , (4.46) 

which yields, 

pl +1 =D 2 L d (ql,ql), p° k = -D^&ql), q\ = q° k+1 . (4.47) 

4.3.4. (-)-Discrete Dirac Structure and Implicit {—^-Discrete Hamiltonian System. Let H d _ ■ TL- — > K be 
given. Then 

dH d .(p° k , q \) = D 2 H d -{p\,ql) dpi + DxHi-bl,^) <kl & T*H- (4.48) 
An implicit (-)-discrete Hamiltonian system (H d _, Aq, X d ) is defined as 

(X k dl dH d _) e D d A ~ <=► («?°, g° +1 ) e Aq and Sl b d _(X$) e A^_ . (4.49) 

Now 

dH d _{ql,pl) = (p k ,ql,D 1 H d .,D 2 H d _), (4.50) 

and 

= (p2 >g g +11 - P ° k+1 ). (4.51) 

So we obtain the set of equations 

ql = -D 1 H d _(jpl,ql), p° k+1 + D 2 H d -(jP k) <£) € A Q (^), ^ = (gg, gg +1 ) £ A Q , (4.52) 

which we shall call the implicit (-)-discrete Hamilton's equations. 

If Aq = Q x Q, the discrete Dirac structure is given by (|4.46|) . and the implicit (— )-discrete Hamilton's 
equations have the form 

<& = -D 1 H d -(j%,ql), pl +1 = -D 2 H^(p° k ,ql), q\ = q° k+1 . (4.53) 

Note that this is essentially the (— )-discrete Hamilton's equations in Lall and West [20| . 

5. Hamilton-Pontryagin Principle and Implicit Lagrangian Systems 

We will first review the continuous Hamilton-Pontryagin principle in coordinates and without constraints, 
to motivate the construction of the discrete Hamilton-Pontryagin principle. A variational derivation of the 
implicit Euler-Lagrange equations can be obtained by considering the augmented variational principle given 
by, 

5 f[L(q,v)-p(v-q)}=0, (5.1) 

where wc impose the second-order curve condition, v = q using Lagrange multipliers p. This variational 
principle is also referred to as the Hamilton-Pontryagin Principle, and is a variational principle on the 
Pontryagin bundle TQ ®T*Q. 

In a local trivialization Q is represented by an open set U in a linear space E, so the Pontryagin bundle 
is represented by (U x E)®(U x E*) = U x E x E* , with local coordinates (q, v,p). If we consider q, v, and 
p as independent variables, we have that, 



6 [L(q,v)-p(v-q)]dt 



Sq + ( — p I Sv — (v — q)5p + pSq 



dq \ dv 



dt 



dt (5.2) 



where we used integration by parts, and the fact that the variation Sq vanishes at the endpoints. This yields 
the implicit Euler-Lagrange equations, 

. dL dL . 

V = q - (5 - 3) 
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5.1. Variational Principle in Phase Space. By starting with the Hamilton-Pontryagin Principle, and 
considering the necessary stationarity conditions in different orders, we will obtain the variational principle 
on phase space, and the usual variational principle on TQ. 

By taking variations with respect to v, we obtain the relation, 



— (q,v) -p = 0. 
ov 

We introduce the Hamiltonian, H : T*Q — > R, defined to be, 

H (q, p) = pv- L(q, v)\ p=dL /dv(q v) . 
From the definition of the Hamiltonian, we can express the augmented variational principle as, 

5 f [pq-H(q,p)] = 0, 



(5.4) 
(5.5) 
(5.6) 

which is the variational principle in phase space. 

By first taking variations of the augmented variational principle with respect to p, we obtain the usual 
variational principle on TQ, 

(5.7) 



5 J L(q,q)=0. 

5.2. Hamilton's Equations. By considering Hamilton's principle in phase space, we obtain, 

= 6 J [pq-H(q,p)]dt 

r)H 1 

alt 



1 



.. dH . dH / 
qdp + pdq - -77-dq - -^-Op 
oq op 

- p --dq-) Sq+ { q --dpr) Sp 



(It. 



By the fundamental theorem of the calculus of variations, this is equivalent to Hamilton's equations, 

. _ dH . _ dH 

^ dp ' P dq 



(5.8) 



(5.9) 



6. Discrete Generalized Legendre Transform 

In this section, we introduce the discrete Lagrangian and Hamiltonian constraints, and develop the dis- 
crete generalized Legendre transform that relates the two discrete constraint submanifolds, by following the 
approach of §2 of [29j and using the symplectic structures Qt'Qxt*q = —d^d± = d\d± on T*Q x T*Q, and 
the maps kq, and fi^i- 

6.1. Discrete Lagrangian constraints. Let Lj, be a discrete Lagrangian on Ag c Q x Q. The symplectic 
manifold (T*Q xT*Q, &t*Qxt*q = -dX d ±) is defined by the quadruple (T*Q x T*Q,Q x Q,itq x ttq, X d ±) 
and the set 

N d := {(z, z 1 ) £ T*Q x T*Q \ (ttq x ttq) ((z, z 1 )) E A d Q , 

X d+ ((w,^ 1 )) = (dL d ({ttq x ttq) ((z,^ 1 ))) ,T{tt q x ttq) ((w^ 1 ))) , 
V(^V) g r ( , iB i) (T*Q x T*Q) s.t. T^i-KQ x ttq) (K^ 1 )) G T (ff£jX ^ )((2 , 2 i )) A d Q } (6.1) 

is a Lagrangian submanifold of (T*Q x T*Q, VLt'QxT'Q = —d\ d ±) of dimension i dim(T*Q x T*Q), where 
the submanifold 

A Q - (ttq xi q )(JV,)cQxQ (6.2) 

is the discrete constraint distribution on Q x Q called the discrete Lagrangian constraint. Then, the discrete 
Lagrangian L d is a generating function of Nd, since N d C T*Q x T*Q is the graph of (hiQ)^ 1 (dL d ). 
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6.2. Discrete Hamiltonian constraints. Let H d ± be a (±)-discrete Hamiltonian on P d ± C 7i±. The 
symplectic manifold {T*QxT*Q, £It*qxt*q = dxd±) is defined by the quadruple (T*QxT*Q, 7i±, t~h ± , Xd±) 
and the set 

N d± := {(z, z 1 ) e T*Q x T*Q \ r H± ({z, z 1 )) € P d± , 

Xd± {(w^w 1 )) = (dH d± (r n± ((z iZ 1 ))),Tt h± ((w^w 1 ))) , 

V( W V) G T {z . zl) (T*Q x T*Q) s.t. T (z , z1) t h± ((w^w 1 )) £ T Tn±{iz , zl)) P d± } (6.3) 

is a Lagrangian submanifold of (T*Q x T*Q,£It*qxt»q = d\ d ±) of dimension idim(T*<5 x T*Q), where 
the submanifold 

P d ± = (ttq x n Q )(N d± ) C H± (6.4) 

is the discrete constraint momentum space called a discrete Hamiltonian constraint. Then, the discrete 
Hamiltonian H d ± is a generating function of N d ±, since N d ± C T*Q x T*Q is the graph of 1 . 

6.3. Symplectomorphism and the discrete momentum function. Consider the identity map, ip : 
{Pi = T*Q x T*Q, fit = -dX d± ) -> (P 2 = T*Q x T*Q, fl 2 = d Xd ±). Since Pi = P 2 and ft 2 = fii, it follows 
that (p*&,2 = fii, and </? is a symplectomorphism. 

The graph of the symplectomorphism ip is a submanifold of Pi x P2, which is denoted by 

r(^) C P x P 2 . (6.5) 

Let 7 V : T(ip) — > Pi x P 2 be the inclusion and 7ri : Pi x P 2 — > Pj be the canonical projection. As in (|2.1| . we 
define 

fyPl XP 2 = T*Ol - 7T2^ 2 

= 7rJ s (-dA d± )- 7 r 2 *(dx d± ). (6.6) 

Since ^ is symplectic, by Proposition 12. 1[ we have that i*Slp lX p 2 = 0. We can write £lp lX p 2 = —dQ d ±, 
where Q d ± = © Xd± = 7T*Ad± + 7r 2 Xd±- Also, r(y) is a maximally isotropic submanifold with half the 
dimension of P 1 x P 2 = (T*Q x T*Q) x (T*Q x T*Q). 

Given the diagonal map * : {T*Q x T*Q) -> (T*Q x T*Q) x (T*Q x T*Q), we have the one-forms **6 d± 
on (T*Q x T*Q), given by 

**6d+ = **(Ad + © Xd+) = Ad+ + <p*Xd+ = (-podqo +pidqi) + (podqo + qidpi) 

= pidqi + qxdpi = d{piqi) = d(G d + ° pf£, Q y), (6.7) 

and 

= **(A d - © Xd-) = ^d- + W*Xd- = (-Podqo +pidqx) + (-pidqi - q dp ) 

= -podqo - qodpo = d(-poqo) = d{G d - o p^. Q y), (6.8) 

where 

Gd+((qo,Qi)®(Qo,Pi))=PiQi, (6.9) 
Gd- ((qo, 9i) © (po, 9i)) = -Po<7o, (6-10) 
which we refer to as the (±)- discrete momentum functions. 

6.4. Discrete Generalized Legendre Transforms. As we saw in (|6.1[) and (|6.3p . the discrete constraint 
manifold on T*Q x T*Q can be realized as graphs of one- forms onQxQ and H±. These various represen- 
tations are related by the discrete generalized Legendre transform, which is a procedure for constructing the 
submanifold !C d ± of the discrete Pontryagin bundle (Q x Q) © Ti± from a submanifold N d of T*Q x T*Q 
associated with (T*Q x T*Q,Q x Q,itq x 7Tq, A(j±) and with a discrete Lagrangian on Aq cQxQ, as in 
(|6.1|) . We define the submanifold K. d ± to be the image of N d under the map ((ttq x ttq) x t« ± ) o 'I', which 
is given by 

fCd± = ((ttq x ttq) xt h± )o y(N d ) C (Q x Q) x ft±, (6.11) 
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which is the graph of the discrete Legendre transforms ¥L d ± ■ Q x Q — > TL± with respect to the discrete 
constraint distribution Aq C Q x Q. Define the discrete generalized energy E d ± on (Q x Q) © 7i-|- in terms 

of the (±)-discrete momentum functions, and pr x , by 

E d ± = G d ±-L d opr( ± , (6.12) 

or more explicitly, 

Ed+ {(qo,qi) © (<?o,Pi)) =Piqi ~ L d {q ,qi), (6.13) 

£d- ((<?o,gi) © (po,9i)) = -paqo - L d (q ,qi). (6-14) 

These expressions agree with the (±)-discrete Hamiltonians H d ± : TL± — > R, once we appropriately impose 
the discrete Legendre transforms FL d ±. More precisely, the following diagram commutes: 



l0xQffiF ^: (Q x Q)®H± 



n± — 



(6.15) 



The maps 1q x q © ¥L d+ and 1q x q © FL d _ are given by {q ,qi) i-» (<7o, <7i, D 2 L d (q a , qi)) and (90,91) ^ 
(?o,?i, — DiL d (q 0l (71)), respectively Then, the submanifold /Q± may be given by 

= {(qo,qi,pi) G (Q x Q)®H+ I V(<j ,Pi) G W+, (qo,qi) £ Aq is a stationary point of E d +) (6.16) 

= {(?o,gi,Pi) G (Q x Q) ©W+ I (q , gi ) g A Q ,p x = £> 2 L d (g , ?1 )} , (6.17) 

and 

K-d- = {(qo,qi,Po) G (Q x Q)®H- I V(p ,gi) G (<7o,9i) G A Q is a stationary point of E d _} (6.18) 

= {(?o,?i,Po) e(QxQ)®H- I (g , <7i) G Aq, Po = -£>i£d(9o,?i)} • (6-19) 

7. Discrete Hamilton-Pontryagin Principle and Implicit Discrete Lagrangian Systems 

We first relax the discrete second-order curve condition, and consider the augmented discrete variational 
principle given by, 

5J2 [LM,ql)-p k+1 (ql - ql +l )\ = 0, (7.1) 

where we impose the second-order curve condition, q], = q® +1 using Lagrange multipliers Pk+i, and keep 
the endpoints q{) an< 4 q% fixed. We will refer to this discrete variational principle as the Discrete Hamilton- 
Pontryagin Principle. 

If we consider q^, q\ and pk as independent variables, with the condition that Sq® = and Sq^_ 1 = 0, we 
obtain, 

[L d {qlql)-p k+1 (q\-ql + l)\ = £ 9k) + pM 

-[Qi - <&+i]Spk+i + [D 2 L d (qlqD - p k+1 ]5q\} (7.2) 

from which we obtain the implicit discrete Euler- Lagrange equations, 

Pk = -D 1 L d (&Qk)> Pk+i=D 2 L d {ql,ql), ql = q° k+v (7.3) 

7.1. Intrinsic Discrete Hamilton-Pontryagin Principle. We now formulate the discrete Hamilton- 
Pontryagin principle intrinsically, by obtaining the intrinsic expression for the discrete Hamilton-Pontryagin 
sum on ((Q x Q) © T*Q) x ((Q x Q) © T*Q). Prior to doing this, we introduce the natural projection maps 
T{QxQ)®n ± ■ ({Q x Q) © T*Q) x ((Q x Q) © T*Q) -» ((Q x Q) © W ± ), which arc given by 



\QxQ)®H + ■ ((<& ql) &> {qlpk)Aq°k+i,ql+i) ® (q°k+i,Pk+i)) » («g,9k)e(gg,Pk+i)> (7.4) 
•(QxQ)aW- : ((<Z°-i,9Li)©(9fe-i,Pfc),(9fe,9fc)©(9fe,Pfc+i)) ^(^,^)©(p fe ,^), (7.5) 
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and pr£ Q ■ (Q * Q) © T*Q -> T*Q given by 

pr d T + Q : (qlql) © (ql,p k ) ^ (&,p k ), (7.6) 

p4~q ■ (ik-uQk-i) © (?fe-i>Pfe) » (<&-i,Pk)- (7.7) 

Recall the projections p^'Qy ■ T*Q x T*Q — > (Q x Q) © W± given by 

P(t.q)» : (8fc+i.Pfc+i)) ^ (<M+i) © (q°k,Pk+i), (7.8) 

P'iT'Q) 2 '■ ((lk-^Pk),(qLPk+x)) >-> (?jfc-l>9jfc) © (Pfc,9*)- (7-9) 

The projections p^.qy andpr^qXpr^Q can be combined to yield a second projection from ((Q x Q) ©T*Q)x 
((Q x <2) © to (Q x Q) © H± that is distinct from T(q x q- )(bh± , namely, 

P d (T* Q y ° (p4 + *q x Wt + ,q) ■ © (llPk), © (3fc+i,Pfc+i)) (3fc,3fc+i) © fak.Pk+i), 

(7.10) 

P(t.q)» ° (P r T*Q x p4~q) : ((?fc-i>?Li) © {ll-i,Pk), (q° k ,Qk) © (9*,Pfc+i)) ^ (3fc_i,«k) © (Pk,qD- (7-11) 
Now, recall the definitions of the (±)-discrete momentum functions from (|6.9|) and (|6.10l) . 

Gd+((qo,qi)®(qo,Pi)) =pm, (7.12) 
Gd- ((?o, ffi) © (po, ?x)) = -pogo, (7-13) 
and the (i)-discrete generalized energies from (|6.13[) and (|6.14|) . 

Ed+ ((?o,9i) © (go, Pi)) =Pl9i - L d {q ,q 1 ), (7.14) 

-Ed- ((go, gi) © (po, gi)) = -pogo - £d(go, gi)- (7.15) 

We can now intrinsically characterize the (±)-discrete Hamilton-Pontryagin sums in terms of the above 
quantities. By direct computation, one can verify that the (+)-discrete Hamilton-Pontryagin sum can be 
expressed as a functional on 

(4,4+i) = ((<&«£) © (qlpk),(q k + i,ql + i) © (92+i,p*+i)) e ((Q x Q) ©t*q) x «q x q) ©t*q) , 

(7.16) 

which is given by 

^L d {ql,q\) -p k+1 (q\ - q° k+1 ) = J^Pk+iQk+i ~ {Pk+lQk ~ L d(q° k ,ql)) 

= H (°d+ ° P^.Q)a ° (P r T^Q X P^q) ^ E d+ ° ^(QxQ)©H + ) (x+, X+ +1 ). 

(7-17) 

Similarly, the (— )-discretc Hamilton-Pontryagin sum can be expressed as a functional on 

(Sfc-i,**) = ((«k-i,«iLi) © («fc-i.j»fc), (g°,gfe) © (qlPk+i)) e ((Q x Q) © t*q) x ((q x q) © t*q) , 

(7.18) 

which is given by 

Y.^iqUl) -Pk{qU -q° k ) = E-p^-i - {-Pkql-L d {qUD) 

= ( Gd - ° P d (T*Q) 2 (P r T~Q X P r T«q) - E d- 7-(QxQ)eW-) (^fc-U^fc)- 

(7.19) 

The intrinsic implicit (±) -discrete Euler- Lagrange equations are given by 

(p4^q x p4^q) Xd± = rf QxQ)@H± dE d ±. (7-20) 

Since we have previously shown that the discrete Hamilton-Pontryagin principle yields the implicit discrete 
Eulcr-Lagrangc equations, it remains to show that 1|7.20[) in coordinates recover (|7.3p . We compute the 
one-forms on the left hand side of each equation, 

c^Pfc+igfc+i = [pk+i d ql+i + gfe+i^Pfc+i] 
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= [Pkdq° k + q° k+1 d Pk+1 ] + ^ [p k+1 dq° k+1 - Pkdql] 

= X! [Pk<tek + Qk+i d Pk+i] + PNdq° N - p dq°, (7.21) 

and 

= [-Pk+idql - ql-idpk] + Y [Pk+i d ll -Pkdql-i] 

= Y ['Pk+idql - ql-i d Pk] +PNdq 1 N _ 1 -p dq 1 _ 1 . (7.22) 

Since wc only evaluate these one-forms on discrete curves with fixed endpoints, the boundary terms vanish, 
and it is sufficient to show that the expressions in the brackets of (|7.2ip and (|7. 22[) agree with the left hand 
side of (|7.20p . By direct computation, we obtain 

(prp, Q x prp,^j X d+{x+,x+ +1 ) = (xd+ ° (p4+ q x p4+ q )) (x+,x+ +1 ) 

= Xd+ [pr d T t Q i:xi),pr d T t Q {xl +1 )^ 

= Xd+{(ql,Pk),{ql+i,Pk+ij) 

= p k dq k +q° k+1 dp k+1 , (7.23) 

and 

(pr^g x p4~q) Xd-^fe-i^fc) = (xd- ° (p4«q x P r T~o)) ( x k-i> x k) 

= Xd- (prT^Q{Xk-i),PrT^Q( x k)) 
= Xd- ((ql^nPhjAq^Pk+i)) 

= -Pk+idql - ql^xdpk. (7.24) 

Computing the right hand side of (|7.20[) yields 

( T (QxQ)®H + dE d+) ( x k' X k+l) = <%"Wx(Jp + )(iii) 

= d{E d+ {(q k ,ql)®(q° k ,p k+1 ))) 
= d(p k+1 ql -L d {ql,q\)) 

= Pk+idql + qldpk+i - D 1 L d {ql,q\)dql - D 2 L d {ql, q\)dql 

= -D.LdiqlqDdql + (p k+1 - D 2 L d (q° k , £)) dq\ + q\dp k+u (7.25) 

and 

( T *Q*Q)®n- dE d-) i x k-i' x k) = d ( E d- °r { Q xQ)(BH _) (x k _ llXk ) 

= d(E d _((qlql)(B(p k ,ql))) 
= d(-p k ql-L d (qlql)) 

= -p k dql - q a k dp k - D 1 L d {ql,ql)dql - D 2 L d (q° k , q\)dql 
= (-p k - D 1 L d {ql,qD) dq° k - D 2 L d (q° k , q\)dql - q° k dp k . (7.26) 
Equating (j7T2U|) with (|7T2"5)) . and (|7T2^|) with ^TW) yield 

Pkdql + q° k+1 d Pk+1 = -D x L d {ql, q\)dql + ( Pk+1 - D 2 L d {q° k , £)) dq\ + qld Pk+1 , (7.27) 

and 

-Pk+idql - ql_ x dp k = (-p k - DiL d (& q l k )) dq° k - D 2 L d {q° k , qftdql - q° k dp k , (7.28) 
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which both recover the implicit discrete Euler-Lagrange equations, 

p k+1 =D 2 L d (q° k ,ql), p k = -D^&ql), q\ = q° k+1 . (7.29) 

8. Discrete Lagrange-d'Alembert-Pontryagin principle 
The (+)-discrete Lagrange-d'Alembert-Pontryagin principle is given by 

5Y,[L d {qUl)-Pk + M-ql +1 )] 

= Y,{[ D ^(llql)+Pk] 5q° k + [D 2 L d (q° k ,ql)-p k+1 ] 5q\ + [ql - q° k+1 ] Sp k+1 } = 0, (8.1) 

for fixed endpoints <?q and q^, and variations (Sq k ,Sql,Sp k ) of (q k ,ql,p k ) £ (Q x Q) © T*Q such that 
5q k £ Aq(^), 8q\ £ Aq(^), and the discrete constraint (q k ,q k ) £ Aq. 

Similarly, the (— )-discrete Lagrange-dAlembert-Pontryagin principle is given by 

6^[LM,9k)-Pk(<&-i-<&)] 

= Y,{[ D i Ld (ti,<ll)+Pk] Sq° k + [D 2 L d (q° k ,ql)-p k+1 ] 5q\ + [ql_ t - q° k ] Sp k } = 0, (8.2) 

for fixed endpoints and q l N _i and variations {5q^,5q\,5p k ) of {q k ) ,q\,p k ) £ (Q x Q) © H- such that 
Sq k £ Aq(^), Sql £ Aq(ql), and the discrete constraint (q k ,q k ) £ Aq. 

Proposition 8.1. The (±)-discrete Lagrange-d'Alembert-Pontryagin principles for the discrete curve {a^ } , 
£ {Q x Q) © T*Q, yield discrete equations that are given in coordinates by 

D 1 L d (4,q 1 k )+p k eA^(4), D 2 L d (q° k ,ql)-p k+1 eA° Q (ql), ql = q° k+1 , (q° k ,ql) £ A d Q . (8.3) 

These equations are equivalent to the coordinate expressions for implicit discrete Lagrangian systems. 

8.1. Constraint distributions on {(QxQ)®T*Q) x ((Q x Q) © T*Q). We will introduce constraint 
distributions for the intrinsic form of the discrete Lagrange-dAlembert-Pontryagin principle. 

Consider the discrete constraint distribution Aq C Q x Q. Recall from Section T4.2.3I that this induces 
a continuous constraint distribution Aq C TQ. In addition, Aq also induces the submanifolds JC d ± C 
(Q x Q) © U± given by (flTTT)) and (flTTT)]) . which were defined in Section [Ql as 

K d+ = {(q Q ,qi,Pi) £(QxQ)®H+\ (q ,qi) £ A d Q , Pl = D 2 L d (qo, qi )} , (8.4) 

and 

K d - = {(q ,q u p ) € (Q X Q) ®H- \ (q , qi ) £ A Q ,p = -DxL d {q , . (8.5) 

Recall also the projection map T ( Q xQ)en± : ((Q x Q) © T*Q) x ((Q x Q) © T*Q) -> ((Q x Q) © W±), which 
induces a discrete constraint distribution A d ± on (Q x Q) © T*Q by 

Ai± := ( T (QxQ)ew±) (£d±), (8-6) 

that encode (q k ,q k ) £ Aq, (<7fc + i , 9a; + i) £ Aq, and the appropriate (±)-discrete momentum constraints. 
We introduce the natural projection prg : (Q x Q) © T*Q — > Q, given by 

pr^:(«fc.?*)©(flg,l'iN-i)'- > 9fc ) (8-7) 
from which we obtain the discrete constraint distribution B d C ((Q x Q) © T*Q) x ((Q x Q) © T*Q) by 

Bd^Cpr^xj^)" 1 ^), (8.8) 

which encodes the constraint (<7°,9fc +1 ) G Aq. 

Let Cd± be the intersection of the discrete constraint distribution B d with A d ± , which is given by 

C d ± := B d n A d± C ((Q x Q) © T*Q) x ((Q x Q) © T*Q) . (8.9) 

We introduce another pair of natural projections 

P r (TOV : T (((Q x 0) © T *Q) x ((Q x Q) © T *Q)) -^TQxTQxTQx TQ, (8.10) 
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given by 

^ {(qlS q l),( q l,S q l),(q k+1 ,6q k+1 ),(ql +v S q l +1 )) , (8.11) 

-> {(ql6ql),(ql6ql),(q k+1 ,5q° k+1 ),(ql +1 ,6q k+1 )), (8.12) 
which lift Aq x Aq x Aq x Aq to continuous distributions T± on ((Q xQJffi T*Q) x ((Q x Q) © T*Q) 

JF ± := (prf Q)4 ) _1 (A Q x A q x A q x A q ) c T (((Q x Q) © T*Q) x ((Q x Q) © T*Q)) , (8.13) 

that encode the constraints 5q° k G A Q (gg), o"g^ G Aq(o|), 8q k+1 G A Q (gg +1 ), and 5q\ +1 G Aq(^ +1 ). 

Let £?± be continuous distributions on ((Q xQ)ffl T*Q) x ((Q xQJffi T*Q) that are obtained by restrict- 
ing T± to C d± , 

Q± ■= F± n rC d± c T (((Q x Q) © T*Q) x ((Q x Q) © T*Q)) . (8.14) 

As we will see, these constraint distributions arise in the intrinsic formulation of the discrete Lagrange- 
d'Alembert-Pontryagin principle. 

8.2. Intrinsic discrete Lagrange d'Alembert Pontryagin principle. The (+)-discrete Lagrange- 
dAlcmbcrt-Pontryagin principle for a discrete curve {x k } on (Q x Q) ®T*Q with fixed endpoints is given 

by 

S J2 { Gd + ° P C {T'Q)i ° ^ r T*Q X P r ftq) - E d+ ° T(QxQ)®-H + ) (»fc . ®fc+l) 

= X! [(p*rt Q x JWt+q) Xd+ - T* QxQ)em+ dE d+ ^ (x+,x+ +1 ) ■ (w+ k+1 ), (8.15) 

which holds for all w+ k+l G 5+ «,a£ +1 ) C T (((Q x Q) © T*Q) x ((Q x Q) © T*Q)). 

Similarly, the (— )-discrcte Lagrangc-d'Alembert-Pontryagin principle for a discrete curve {x k } on (Q x 
Q) ®T*Q with fixed endpoints is given by 

{ Gd - ° P d (T*Q)i ° (P r T*Q X P4«Q) - O T(QxQ)e«-) (^fc-l^fc ) 

= [ ^P r< T*Q X p4«q) Xd- - T (QxQ)®n- dE d-] (n>k-l> X k) ' KT-l,fc)> ( 8 ' 16 ) 

which holds for all w k _ lk G Q- (x^_ v x^) CT(((Q x Q) ©T*Q) x ((Q x Q) © T*Q)). 

Proposition 8.2. TTie {±)-discrete Lagrange-d'Alembert-Pontryagin principles are equivalent to the equa- 
tions, 

(prp-Q x P r ft Q ) Xd+ (xt,x+ +1 ) • (u>£ fc+1 ) = r ( * QxQ)ffiW+ • (w+ fc+1 ) , (8.17) 
/or all w+ k+1 G jZfc+i)* and 

(p4"q x p4-q) Xd- {x k -i,x k ) ■ (u>k-i,k) = T {Q-KQ)<sn- dE d- {x k -i,x k ) ■ (W-iJ » ( 8 - 18 ) 
for all w k l k G (/_ (^fe-ii^fe )> respectively. 

8.3. Intrinsic form of implicit discrete Lagrangian systems. We first introduce (±)-discrete vector 
fields on M, denoted by Xd±{M), to be maps Xd± : M — > M x M, such that 7r^ o Ad+ = 1m and 

o = ljvf, where 7r]^ : M X M —> M; (mo, mi) i— > m^. 
Let Ac(± be (±)-discrete vector fields on P±, where P± = FL^±(Aq), and let Xd± be (±)-discrete vector 

fields whose images are in Ad± = { T (QxQ)®n±) (fcd±)- Consider the discrete integral curves {^}on T*Q 
of X d ±, where z k = (q k ,p k ), i.e., {z k ,z£ +1 ) = X d +{z k ) and (z k _ 1 ,z k ) = X d -(z k ). We consider discrete 
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curves {x^} on (Q x Q) ®T*Q, where x k = {q k ,ql) © (q k ,Pk), and x k = {q k ,q k ) © (q k ,Pk+i), such that the 
second component is given by the curves {z^} on T*Q. Explicitly, we have 

4 = (9fc.9fc)©(9fc.Pfc). ( 8 - 19 ) 
*fc = (9fc.9* + i) © ( 8 - 2 °) 

Using the natural projections pr^q : (Q x Q) © T*Q — > T*<3, which we previously introduced, we impose 
pointwise constraints on the first component, 

(pr d T ± Q x pr d T %) (X d± (xt)) = X d± (p4%(xt)) , (8.21) 

where X,i+(xq) = (xq,Xi), X d —(x~^) = (xZi,Xq), but the curve is otherwise arbitrary. Since we restrict 
the image of X d ± to A d ± , we have the property, 

(p4% x p4± Q ) (X d± {xt)) = X d± (p4%(4)) ■ (8-22) 

Notice that while (|8.21j) give pointwise conditions at k = 0, when we restrict our lifted (±)-discrete vector 
fields to Ad±, we obtain (|B"22|l that holds globally. If {xf }, x\ = {q\,q\) © (q k ,Pk) 6 (<9 x Q) © T*Q, 
• T fc = (#fc-i> 9fc-i) ® (ql-nPk) £ (<9 x <9) © r*<5, are discrete integral curves of X<j±, then it follows that 

X d+ (xt) = {(q k ,q 1 k )®(q k ,p k ),(q k+1 ,q 1 k+1 )®(q° k+1 ,Pk+l)), (8.23) 
= ((<jg_ ls <£_i) © (^_i,p fc )» (9fc»«fc) © (ff£.Pk+i)) • (8-24) 

Proposition 8.3. Let {x k } be discrete integral curves of the (±)-discrete vector fields X d ± on (Q xQ)(BT*Q 
that are naturally induced from (±)-discrete vector fields X d ± on T*Q. If {x k } are discrete solution curves 
of the (±)-discrete Lagrange-d'Alembert-Pontryagin principles, then they satisfy 

(pr d + Q x P 4+ Q y X d+ (X d+ (x+)) ■ (w+ ik+l ) = T( QxQ)en+ dE d+ (X d+ (.x+)) ■ (w+ k+1 ) , (8.25) 
for all w+ k+1 e G+ (Xd+ix^yj, and 

(p4~Q x P r T~ Q ) Xd- (Xd-ixl )) • (Wk-l,k) = T (QxQ)®H- dE d- (Xd-(x k )) ■ («V-l,fc) > ( 8 ' 26 ) 

for all w k _ hk E Q- (JjC d -(x k ty . 

Proof. Since {xf} are discrete integral curves of X d ±, respectively, we obtain (|8.23[) and (|8.24[) . Substituting 
the respective equations into (|8. 1T[) and (|8.18|) yields the desired result. □ 

We may now summarize our results in the following theorem. 

Theorem 8.4. Let L d be a discrete Lagrangian on Q x Q and Aq be a discrete constraint distribution on 
Q. Let X d ± be (zc)- discrete vector fields on P± = ¥L d ±(Aq), such that (L d , Aq, X d ±) are implicit (±)- 
discrete Lagrangian systems. Denote discrete curves on (Q x Q) ®T*Q by {x^}, where x k = (q k , q k ,p^) 6 
(Q x Q) ®T*Q, x k = (q k , q k ,P k+1 ) £ (Q x Q) ©T*Q. Then, the following are equivalent: 

(i) {x k } are discrete solution curves of the implicit (±)-discrete Lagrangian systems (L d , Aq, X d ±). 

(ii) {x k } satisfy the (±)-discrete Lagrange-d'Alembert-Pontryagin principles. 

(iii) {x k } are the discrete integral curves of (±)-discrete vector fields X d ± on (Q x Q) © 7i± that are 
naturally induced from X d ± . 
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9. Discrete Hamilton's phase space principle 
By taking variations with respect to q\ and q® in the discrete Hamilton-Pontryagin principle, we obtain 

D 2 L d (q° k ,ql)-p k+1 =0, (9.1) 
D 1 L d (q° k ,q 1 k )+p k = 0, (9.2) 

respectively. These equations define the discrete fiber derivatives, IFX^ : Q x Q — > T*Q, which are the 
discrete analogues of the Legendre transform, and are given by 

FL+ : (qlql) ^ (ql D 2 L d (q° k , «£)) , (9.3) 

¥L- : ( q ° k , q\) ^ (q° k , -DiL d {&,£)) . (9.4) 

When the discrete fiber derivatives are invertible, we may construct the lift : H± — > /Cd±, where we let 
Ag = QxQ. Since £d± C (Q x Q)©W±, we may, by a slight abuse of notation, let : H± — > (Q x Q)©W±. 

We introduce the (±)- discrete Hamiltonian, H d ± ■ H.± — > K, which we define in terms of the discrete 
generalized energy E d ± : (Q x Q) ffi W± — > R to be 

H d± =E d± oa d ±, (9.5) 

which agrees with the usual definition of the (±)-discrete Hamiltonians given by 

H d+ {q° k ,p k+1 ) = Pk+iql - Ld{q° k ,qi)\ pk+i=D2Ld{qUl} , (9.6) 

H d -(pk,q° k+1 ) = ~p k q ( k ~L d (q° k ,ql)\ pk= _ DiLd{q „ q , y (9.7) 

or equivalently, 

Hd+{qk,Pk+i) = p k +iq k +i - L d{qk,qk+i)\ Pk+l=D2Ld ( qk ^ qk+l ) , (9.8) 

H d ^{p k ,q k+1 ) = -p k q k - L d (q k ,q k+1 )\ pk= _ DiLd(qk qk+i) . (9.9) 

From the definition of the discrete Hamiltonian, we can express the discrete variational principle as 

Sj2\Pk+iqk+i - H d+ (q k ,p k+1 )} = 0, (9.10) 

5J2[-Pkqk-H d _(p k ,q k+1 )]=0, (9.11) 

where we keep the endpoints go and q^ fixed. This gives the (±)- discrete Hamilton's principles in phase 
space. 

9.1. Discrete Hamilton's equations. By considering the (±)-discrete Hamilton's principle in phase space, 
we obtain, 

= 8^2\p k+1 q k+ i - H d+ (q k ,p k+1 )} 
= y^{[gfc+i - D 2 H d+ {q k ,p k+ i)} Sp k+ i + [p k - DiH d+ (q k ,p k+ ij\ 5q k } , (9-12) 

and 

= S^2[-p k q k - H d -(p k ,q k+ i)] 
= / ]{[-qk - DiH d -(p k ,q k+ i)] 5p k + [-Pk+i - D 2 H d -(p k ,q k+1 )} Sq k+1 } , (9.13) 
From which we obtain the (+)- discrete Hamilton's equations, 

q k+ i = D 2 H d+ (q k ,p k+1 ), p k = D-yH d+ (q k ,p k+1 ), (9.14) 

and the (-)-discrete Hamilton's equations 

q k = -D 1 H d _(p kl q k+1 ), p k+1 = -D 2 H d _(p k , q k+ i). (9.15) 
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9.2. Intrinsic discrete Hamilton's phase space principle. We recall the two pairs of maps, P?y.Q)2 : 
T*Q x T*Q -> (Q x Q) © W ± , and r w± : T*Q x T*Q -» H±. The intrinsic (+)-discrete Poincare-Cartan 
sum is given by 

^2 \Pk+lQk+i - #d+(<7fc,Pfc+i)] = ^2 [ Gd + ° P d (T*Q) 2 ~ H d+°T U+ (zk,Zk+i), (9.16) 
and the intrinsic (— )-discrete Poincare-Cartan sum is given by 

~ H d-{Pk,qk+i)} = ^2 [ Gd - °P d (T'QY- ~ Hd - OT n- {Zk,Zk+l), (9.17) 

where z k = (qk,Pk)- We view both of these sums as functionals on (z k , Zfe+l) G T*Q x T*Q. 
The intrinsic (±)-discrete Hamilton's equations are given by 

Xd±=r^ ± dH d± . (9.18) 
Computing the exterior derivative of the discrete Poincare-Cartan sums yield 
[Pk+iQk+i - H d+ (q k ,p k+1 )} 

= ^ [pk+idqk+i + qk+idpk+i - dH d + o Tn + (zk, Zfc+l)] 

= ^2 [Pkdqk + Qk+idpk+i] - ^ T H+dH d +(zk, z k +i) + ^ [p k +idq k+1 - p k dq k ] 

= X! iXd+(z k , z k+1 )} - ^2T^ + dH d+ {z k , Zk+i) + [pNdq N ~ Podq ] 

= X! [*d+ ~ T H+dH d+ ~\ {zk, Zfc+i) + [ptv^w - Podqo] , (9.19) 

and 

t~P fe * ~ H d-(Pk,qk+i)} 

= ^2 [^Pkdqk - qkdpk - dH d - o t h _ (z k , Zfe+i)] 

= ^ [-pfc+idgfe+i - qfcdpfc] - i~Z i _dH d -(z k , z k+1 ) + ^2 [Pk+idqk+i - Pkdqk] 

= ^2 [Xd-(z k , z k+1 )} - (z fc , Zk+i) + [pNdq N - Podqo] 

= ^2 [xd-Tn_dH d ^ (zk, Zfc+i) + [pNdq N - Podqo] ■ (9.20) 

Since we only evaluate these one-forms on discrete curves with fixed endpoints, the boundary terms vanish, 
and we recover (|9.18[) . 

10. Discrete Hamilton-d'Alembert principle in phase space 

We show how the implicit discrete Hamiltonian systems can be derived from a generalization of the 
discrete Hamilton's principle in phase space, which we refer to as the discrete Hamilton-d'Alembert principle 
in phase space. 

The (+)-discrete Hamilton-d'Alembert principle in phase space for a discrete curve {(q k ,p k )} in T*Q is 
given by, 

S ^2lPk+iqk+i - H d+ (q k ,p k+1 )] = 0, (10.1) 
where we require that (q k , q k +\) G Aq, and the variation, when the endpoints are kept fixed, is given by 

3^2[Pk+iqk+i - H d+ (q k ,Pk+i)] 

= X! - D 2H d+ (q kl p k+1 )] Sp k+1 + [p k - D 1 H d+ (q k ,p k+1 )] 5q k } , (10.2) 

where we choose constrained variations 6q k G Aq(q k ). 

The (— )-discretc Hamilton-d'Alembert principle in phase space for a discrete curve {(q k ,p k )} in T*Q is 
given by, 

SY^l-PhOk - H d -(q k ,p k+1 )] = 0, (10.3) 
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where we require that (q k , qk+i) G Aq, and the variation, when the endpoints are kept fixed, is given by 

S^2[-PkQk - H d _(p k ,q k+1 )) 

= ^{[-Qk - D 1 H d ^{p k ,q k+1 )]5p k + [-Pk+i - D 2 H d -(p k ,q k+ i)}6q k+ i} , (10.4) 

where we choose constrained variations 6q k £ Aq ((?&). 

Note that if one starts with the (±)-discrcte Lagrangc-d'Alembcrt-Pontryagin principle, and optimize 
over q\ and q^, respectively, one obtains the (±)-discrctc Hamilton-d'Alembert principle in phase space. 

Proposition 10.1. The {±)-discrete Hamilton-d'Alembert principle in phase space for a discrete curve 
{{lk,Pk)} in T*Q give the implicit (±)-discrete Hamiltonian systems in (|4.32p and (|4.49[) , respectively. 

Proof. From (|10.2p . the (+)-discrete Hamilton-Pontryagin principle is equivalent to 

[q k+ i - D 2 H d+ (q kl p k+ i)} Sp k+ i + [p k - DiH d+ (q k ,p k+ i)] Sq k = 0, (10.5) 

for all 5q k e Aq(%), all Sp k+ i and all (q k ,q k+ i) £ Aq. Thus, we obtain (|4.35[) which are (|4.32p expressed 
in coordinates. Similarly, from (|10.4[) . the (— )-discrcte Hamilton-Pontryagin principle is equivalent to 

[-q k - D 1 H d -{p k ,q k+ i)} Sp k + [—Pk+i - D 2 H d -(pk,qk+i)] Scj k+1 = 0, (10.6) 

for all Sq k+ i £ Aq((7a; + i), all 5p k , and all (q k ,q k+1 ) £ Aq. Thus, we obtain (|4.52|) which are ()4.49[) expressed 
in coordinates. □ 

10.1. Constraint distribution on T*Q x T*Q. As with the intrinsic formulation of the discrete Lagrange- 
d Alcmbcrt-Pontryagin principle, we require a constraint distribution on T*Q x T*Q to formulate the discrete 
Hamilton d'Alembert principle intrinsically. 

Define a discrete constraint distribution on T*Q by 

A d T , Q := (ttq x ttq)- 1 (A d Q ) C T*Q X T*Q, (10.7) 

where 7Tq : T*Q — ► Q. Recall the discrete fiber derivatives, FL^ : Q x Q — > T*Q. Then, we let P± be the 
image of the discrete constraint distribution Aq under the discrete fiber derivatives, i.e., P± = FL^ (Aq) C 
T*Q. Let Ap ± be the restriction of A t „q to P±, which is given by 

A d P± := A^q n (P± x P±) c T*Q x T*Q. (10.8) 

We canonically identify T(Q x Q) with TQ x TQ and, by a slight abuse of notation, write T(ttq x ttq) : 
T (T*Q x T*Q) TQ x TQ. Define the continuous constraint distribution on T*Q x T*Q by 

1 := (T(ttq x ^q))- 1 (Aq x A q ) C T (T*Q x T*Q) . (10.9) 

Then, we let J± be the restriction of I to A P± , given by 

J±:=lnTA d P± cT(T*QxT*Q). (10.10) 

10.2. Intrinsic discrete Hamilton-d'Alembert principle. The (±)-discrete Hamilton-d'Alembert prin- 
ciples for discrete curves {z^} on T*Q with fixed endpoints are given by 

SJ2 [G d ± o pfi» Q)2 - H d± o r H± ] = ]T [ Xd ± - r^ ± dH d± ] (z±, z± +1 ) ■ ( S ± fe+1 ) , (10.11) 

which holds for all s^ fe+1 £ J±(zjjr, zf +1 ) C T (T*Q x T*Q). This yields the following two propositions. 
Proposition 10.2. The (±)-discrete Hamilton-d'Alembert principles are equivalent to the equations, 

Xd± (z k ,z k+1 ) ■ (sfc,fc+i) = Tn ± dH d ± {z^,z^ +1 ) • (s± fe +i) . ( 10 - 12 ) 
for all s± fc+1 £ J±(z±, z± +1 ). 
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Proposition 10.3. Let {z k } be discrete integral curves of the (±)-discrete vector fields X d ± on T*Q. If 
{z^} are discrete solution curves of the (±)-discrete Hamilton-Pontryagin principles, then they satisfy 

Xd+ (X d+ (z+)) • (< fc+1 ) = r^ + dH d+ (X d+ (z+)) ■ (,s+ fe+1 ) , (10.13) 
for all s^ k+1 e J+ (X d+ (z£)), and 

Xd - (X d ^(z~ +1 )) • (sfc ifc+1 J 

for all s k k+1 6 J71 (X d - 1)) > respectively. 

Proof. Since {z k } are discrete integral curves of X d ±, we obtain, X d+ {z k ) = (z k ,z k+1 ) and X d ~(z k+1 ) = 
(Zk' z k+i)> wn i cn together with (|10.12[) yield the desired results. □ 

We may now summarize our results in the following theorem. 

Theorem 10.4. Consider (±)-discrete Hamiltonians on TL± and a given discrete constraint distribution Aq 
on Q. Let X d ± be (±)- discrete vector fields on T*Q, such that (H d ±, Aq, X d ±) are implicit (±)- discrete 
Hamiltonian systems. Let {z k } = {(l^ ,Pk)} be discrete curves onT*Q. Then, the following are equivalent: 

(i) {z k } are discrete solution curves of the implicit (±)-discrete Hamiltonian systems (H d ±, Aq, X d ±). 

(ii) {z k } satisfy the (±)-discrete Hamilton-d'Alembert principles in phase space. 

(iii) {z k } are the discrete integral curves of the (±)-discrete vector fields X d ± on T*Q. 

11. Conclusion 

In this paper, we developed the theoretical foundations of discrete variational mechanics arising from 
the discrete Hamilton-Pontryagin principle, which is a discrete variational principle that provides a unified 
treatment of discrete Lagrangian and Hamiltonian mechanics, as well as the discrete analogue of Dirac 
structures that unifies symplectic and Poisson structures. We demonstrate that just as variational integrators 
preserve a discrete symplectic structure, Hamilton-Pontryagin integrators preserve a discrete Dirac structure. 

Discrete Dirac structures are intimately related to the geometry of Lagrangian submanifolds, which arise 
in the geometric description of the Hamilton- Jacobi equation and the theory of canonical transformations. 
Indeed, our construction of a discrete Dirac structure is given in terms of the maps Kq : T*QxT*Q — > T*(Qx 
Q) and fl d ± : T*Q x T*Q — > T*7i±, that are obtained naturally from the geometric description of generating 
functions of symplectic maps. This yields a discrete geometric formulation of implicit discrete Lagrangian and 
Hamiltonian systems, thereby providing a unified theoretical foundation for developing geometric numerical 
integrators for degenerate, interconnected, and nonholonomic Lagrangian and Hamiltonian systems. 

We introduced an extended discrete variational principle, called the discrete Hamilton-Pontryagin princi- 
ple that is valid for all discrete Lagrangians, and provide an intrinsic description that is valid semi-globally. 
Discrete constraints have been incorporated by the introduction of the discrete Lagrange-d'Alembert- 
Pontryagin principle, and the corresponding discrete Hamiltonian description is given in terms of the discrete 
Hamilton's principle in phase space, and the discrete Hamilton-d'Alembert principle in phase space. These 
discrete variational principles establish a link between discrete variational structures, discrete Dirac struc- 
tures, and implicit discrete Lagrangian and Hamiltonian systems. 

While this paper is motivated by the desire to understand the geometry of Dirac integrators in the context 
of geometric numerical integration, it is worthwhile to note that the discrete Lagrangians and Hamiltonians 
in our theory are simply generating functions of Types 1, 2, 3, and as such, the resulting theory of discrete 
Dirac mechanics provides a general characterization of near- identity Dirac maps. 

Several interesting topics for future work are suggested by the theoretical developments introduced in this 
paper: 

• Discrete Dirac formulation of the discrete Hamilton-Jacobi equation with a view towards de- 
veloping a discrete analogue of the Hamilton-Jacobi theory for nonholonomic systems [l5| . Since 
discrete Dirac structures are related to Lagrangian submanifolds, which in turn describe the geometry 
of the Hamilton-Jacobi equation, it is natural to explore the Dirac description of Hamilton-Jacobi 



= r* H _dH d _ (A d _(z- +1 )) • ( S - fc+1 ) , (10.14) 
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theory, as Dirac structure can incorporate nonholonomic constraints and thereby provide a unified 
treatment of both the classical and nonholonomic Hamilton-Jacobi theory. 

• Galerkin variational Hamiltonian integrators for degenerate systems, by a careful application of the 
generalized discrete Legendre transformation to Galerkin variational Lagrangian integrators. The 
Galerkin approach to constructing symplcctic methods on the Lagrangian side leverages results in 
approximation theory to obtain generalizations of symplectic methods that incorporate adaptive, 
multiscalc, and spectral techniques (Chapter 5 of By embedding the theory of Galerkin varia- 
tional integrators into discrete Dirac mechanics, and considering the Hamiltonian analogue, we will 
obtain a general framework for constructing symplectic methods with prescribed numerical approx- 
imation properties for degenerate Hamiltonian systems, such as point vortices (26|. 

• Discrete reduction theory for discrete Dirac mechanics with symmetry. The Dirac formulation of 
reduction provides a means of unifying symplectic, Poisson, nonholonomic, Lagrangian, and Hamil- 
tonian reduction theory, as well as addressing the issue of reduction by stages. The discrete analogue 
of Dirac reduction will proceed by considering the issue of quotient discrete Dirac structures, and 
constructing a category containing discrete Dirac structures, that is closed under quotients. 

• Discrete multi-Dirac mechanics for Hamiltonian partial differential equations. Dirac generalizations 
of multisymplcctic field theory, and their corresponding discretizations will provide important insights 
into the construction of geometric numerical methods for degenerate field theories, such as the 
Einstein equations of general relativity. 
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